Light curve analysis and evolutionary status of four newly identified short-period eclipsing binaries

We present the physical and orbital parameters of four short-period eclipsing W UMa systems: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ZTFJ000030.44 +391106.9$$\end{document}ZTFJ000030.44+391106.9 (referred to as S1), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ZTFJ000817.08+402532.1$$\end{document}ZTFJ000817.08+402532.1 (referred to as S2), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ZTFJ002158.44+252934.04$$\end{document}ZTFJ002158.44+252934.04 (referred to as S3), and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ZTFJ003357.62+415747.8$$\end{document}ZTFJ003357.62+415747.8 (referred to as S4). The absolute parameters and evolutionary status of these systems are determined, and new times of minima are calculated. Additionally, we present the 3D fill-out configuration for each system. The four Systems exhibit moderate contact W UMa binary with a fill-out factor of 49%, 38%, 28%, and 51%, respectively. Comparing the systems’ periods, we observed a proportional relationship, where shorter periods correspond to lower fill-out factors, and longer periods were associated with higher fill-out factors. Based on the derived surface temperatures and mass ratios of the components, all systems are classified as A-type W UMa binaries. The obtained parameters in addition to a list of previously published data are then utilized to derive an updated Mass-Luminosity relation (M-L) for both A and W-type eclipsing W UMa systems. A comparison with previously published relations reveals that the majority of the EW systems lie between 0.2 and 2 Msun on the M-L diagram. Moreover, we discuss the dynamical evolutionary aspects and evolutionary status of the four components, along with their positions on the Zero Age Main Sequence (ZAMS) and Terminal Age Main Sequence (TAMS).


Data
The photometric observations in this study were obtained from the ZTF database, accessible at https:// www.ztf.calte ch.edu/ ztf-public-relea ses.html, using the gri filters.The collected light curves are generated by using measurements from calibrated single-exposure PSF-fit-derived catalogs.This approach guarantees that there is no contamination of light from nearby systems.Details about the ZTF data reduction can be found at 27 .We conducted a search in the ZTF catalog for variable stars as introduced by 28 .Subsequently, we visually inspected individual light curves of eclipsing binaries and gathered photometric data for systems meeting specific criteria, including: (1) Systems that were conclusively identified as eclipsing binary systems by 27 .(2) Systems that had not been previously analyzed.(3) Systems with orbital periods shorter than 0.3 days, selected because they could be potential contact binaries.(4) Systems for which the visual magnitude and coordinates were suitable for follow-up observations using the Kottamia Astronomical Telescope (a 1.88 m telescope operated by NRIAG).( 5) Systems with light curves exhibiting well-defined photometric variations.After collecting data from the ZTF catalogue, it underwent careful curation, which involved visual inspection to identify and subsequently remove outlier data points.Specific details about the selected objects are available in Table 1.The period of the four systems is determined by 28 .To test and refine the period estimation, we applied the Lomb-Scargle Periodogram method 29,30 on the collected data.
For each of the four systems, we determined the primary and secondary minima timings using the method described by 31 .Our calculations yielded four primary and four secondary minima for each system, as shown in Table 2.We then computed their average values to establish the ephemeris formula.The average primary and secondary minima are provided in Table 3. Table 3 also presents the previously estimated periods alongside those estimated in this study, with discrepancies typically occurring in the fifth decimal place.The phase magnitude diagram for each system is shown in Fig. 1.These values hold significance for future research, enabling the calculation of period changes and the construction of O-C diagrams.Subsequently, we employed the times of primary minima (epoch) ( HJD_Min ) to formulate the ephemeris equations (1, 2, 3, and 4) for each system.
The symbol "E" in the equations denotes the number of integer cycles.

Light curve analysis
The light curve analysis of the four systems was performed using PHOEBE code 32 .To determine the surface temperature of the primary component ( T 1 ) for each system, we followed a two-step process.Firstly, we estimated the color index J-H (e.g., from the 2MASS catalogue) and determined its corresponding temperature using the color-temperature calibration method outlined in 33 .Secondly, we integrated temperature data from Gaia DR3.Finally, we computed the mean value of these temperature estimates.This temperature was fixed through the light curve modeling.Gravity darkening and bolometric albedo exponents for the effective temperature below 7500K (i.e.convective envelopes) were fixed at g 1 = g 2 = 0.32 34 and A 1 = A 2 = 0.5 35 , respectively.Logarithmic limb darkening coefficients (X and Y) were derived from tabulated values using the method of 36 .Due to the lack of spectroscopic observations, q-search method (see, 37 ) was used to estimate the initial value of the mass ratio ( q = M 2 /M 1 ) for each system and then adjusted during the light curve fitting.Figure 2 shows the relation between the residual of input parameters and the mass ratio q (for mass-ratio q ranging from 0.1 to 0.9), where a minimal residual occurs at q= 0.31, 0.13, 0.40 and 0.43 for the systems S1, S2, S3 and S4, respectively which were adopted as the initial mass ratio.In this q scan, the primary star's effective temperature ( T 1 ) and the mass ratio q remained fixed, while the secondary star's effective temperature ( T 2 ), orbital inclination i, the modified potential ( 1 = 2 ), and the passband luminosity of the primary ( L 1 ) were adjusted.According to the light curve shape (EW), we selected the "Overcontact not in thermal contact" mode in PHOEBE.parameters during the fit are T 2 , , i, q and L 1 .The best fit from the PHOEBE code for each system is shown in Fig. (3).The error associated with these parameters is obtained from PHOEBE code at the best fit model.Table 4 lists the integrated (gri) parameters extracted from the best fit of the light curve modeling for each system.Table 2. Times of primary and secondary minima of the four systems.

Star ID Minima I Minima II
The W UMa systems in general are classified into two types based on the temperature and mass of their components.In A-type systems, the hotter and brighter star is more massive than its companion, while in W-type systems, the hotter companion is less massive [38][39][40] .According to the resulting parameters and Following 38 , all the systems were categorized as A-type EW stars.The estimated mass ratio of our sample validates the findings presented by 41 and 42 , indicating that A-type EW stars generally exhibit a lower mass ratio ( q < 0.72 ) than W-type systems.O'Connell effect (asymmetry in the two maxima of the light curve) 43 was detected only in S4 with a hot spot on the surface of the secondary component.5][46][47] ).The surface temperature of the spot is found to be 20% higher than the surrounding photosphere.
We determined the fill-out factor (f) for our systems using the method described by 49 , which is given by the equation f = (� inner − �)/(� inner − � outer ) .In this equation, outer and inner denote the outer and inner criti- cal equipotential surfaces, respectively, while represents the equipotential surface of the star.Our calculations yielded fill-out factor values of 49% ( ±2 ), 38% ( ±1 ), 28% ( ±2 ), and 51% ( ±2 ) for the four systems under study.These results suggest a moderate degree of contact configuration for these systems.For a visual representation of the system's Roche geometry, refer to Fig. (4), which was generated using Binary Maker (BM3) 50 .Notably, we observed a proportional relationship between the fill-out ratio and the orbital period; systems with shorter periods exhibited lower fill-out ratios, and vice versa.Evidence of such correlation is observed when plotting the well-studied W UMa systems reported in 48 (see Fig. 5).This result may lend support to the explanation of a period cut-off observed in this class of eclipsing binaries, as reported by 21 .
In order to trace the evolutionary status for the presented systems, the absolute parameters for each system are calculated.

M-L relation for EW systems
Several attempts have been made to get the mass luminosity relation of eclipsing binaries.These determinations from the early twentieth century see 51 for a historical review. 52constructed a relation for EW type in particular. 51sed a sample of 268 detached and double lined binary stars with accurate T eff and luminosity measurements to improve the M-L relation.They divided the masses into four groups of low mass, intermediate mass, high mass and very high mass stars and each group has its mass-luminosity relation.Recently, 48 made a statistical study on ≈ 700 W UMa eclipsing binaries with published stellar parameters and concluded two equations for the mass luminosity relation for both the primary and secondary from linear fits (see their equations 9 and 10).
In the present work, we propose an update for these M-L relations for the EW UMa systems.Our sample is collected from previously published absolute parameters including 48 and 53 .The selected sample is characterized by being spectroscopically observed or systems in total eclipse.This will help reduce the uncertainty in the estimation of the mass ratio and consequently the rest of the absolute parameters.Table 5 lists the first five rows of the data used in the present work, while the full version of the table can be found at https:// doi.org/ 10. 5281/ zenodo.10209 449.
Our linear fits resulted in the following relations, Figure 6 reveals that the majority of absolute mass and luminosity values for W UMa components are concentrated within the range of 0.1 to 1.5 M sun for the secondary component, whereas for the primary components, they are found to be within the range of 0.6 to 2 M sun .Compared to 51 , our results occupied the low and inter- mediate M-L domains.

Absolute parameters and evolutionary status
To determine the absolute parameters of our systems, four empirically established mass-period relationships (M-PR) were employed to assess the mass of the primary star in W UMa binary systems.The initial M-PR was introduced by 54 , while subsequent relationships were developed by 55 and 56 , and more recently by 48 .The average mass obtained from the four relations is designated as the mass ( M 1 ) of the primary component, subsequently, we derived the mass of the secondary star based on the mass ratio obtained through light curve modeling.The semi-major axis ( a ) in solar radii ( R sun ) was determined using Newton's formulation of Kepler's third law, where (M 1 + M 2 ) = 0.0134 • a 3 /P 2 .With the effective radius ( r i ) obtained from our light curve analysis and the estimated semi-major axis, we calculated the radii of both the primary and secondary stars in solar units as R = r i • a .To estimate the luminosities of the primary ( L 1 ) and secondary ( L 2 ) stars in solar units, we employed www.nature.com/scientificreports/ the Mass-Luminosity (M-L) relation presented in this study.Additionally, we computed the surface gravity ( log(g) ) of the components using Eqs.( 4) and ( 5) introduced by 9 .The parameter uncertainties were determined by considering the error bars associated with the relevant parameters.
Utilizing the parameters listed in Table 6, we represented the evolutionary status of our systems on the Mass-Luminosity (M-L) and Mass-Radius (M-R) diagrams for the Zero Age Main Sequence (ZAMS) and Terminal Age Main Sequence Stars (TAMS).These diagrams were constructed based on the evolutionary tracks provided by 57 with a metallicity (Z) of 0.014.The upper panel of Fig. 7 illustrates the positions of the four systems on the M-R track, while the lower panel of Fig. 7 displays the components on the M-L track.These figures clearly indicate that the primary stars in our systems primarily reside on the main sequence, indicating their relatively less evolved nature.Conversely, the secondary components of the four systems appear above the main sequence, signifying their evolved nature.Notably, the secondary component of S4 exhibits a relatively less evolved behavior compared to its counterparts, while S2 appears to be the most evolved component among them.Furthermore, we examined the dynamical evolution of the systems by studying the orbital angular momentum and total mass of the binary in relation to the orbital period, as investigated by 58 .We employed their Eq.(1) to calculate the orbital angular momentum for our systems.Additionally, we utilized Table 2 to construct the log(J o )-log(P) and log(M)-log(P) diagrams, as shown in Fig. 8.Our findings indicate that the four systems are positioned towards the lower left corners in both diagrams, as expected for W UMa stars.

Conclusion
In this study, we conducted a comprehensive analysis of four short-period eclipsing W UMa systems based on ZTF observations.By deriving absolute parameters such as luminosity, radius, mass, and surface gravity, we gained valuable insights into the physical properties and evolutionary status of these systems.The positioning of the systems on Mass-Luminosity and Mass-Radius diagrams indicated that the primary stars mainly reside on the main sequence, while the secondary components exhibit signs of evolution.Furthermore, the confirmation of the systems' classification as W UMa binaries through dynamical evolution analysis strengthens our understanding of these systems.The detection of the O'Connell effect in one system (S4) allowed us to estimate spot parameters, revealing a hot spot on the surface of the secondary component, which could expand by mass transfer from the  primary component to the secondary one.we also found that, the primary surface temperature as well as mass is higher than the secondary one, indicating A-sub type of W UMa systems.We updated the Mass-Luminosity relations through collecting a well studied sample with spectroscopic and/or totally eclipsed stars.This study also highlighted a relationship between the fill-out factor and the period, however this needs further investigation as well as observations, particularly near the period cut-off (i.e. less than 0.3 d ).The evolutionary status analysis demonstrated a noticeable evolution in the secondary components compared to the primary components in each system.Finally, we recommend conducting spectroscopic follow-up observations of these systems for further investigations and deeper insights into their properties.

Figure 1 .
Figure 1.Phase-magnitude plots of the four systems at gri filters.

Figure 2 .
Figure 2. Relation between the mean residual of input parameters and the mass ratio q of the four systems.The dashed red line refers to the selected initial mass ratio.

Figure 3 .
Figure 3. Best-fitted models for the four systems (S1, S2, S3, and S4) are shown.The solid lines represent the synthetic data extracted from the PHOEBE code, while the colored dots depict the observations.The lower panel displays the residuals for each system.The X and Y axes represent the phase and magnitude, respectively.

Figure 4 .
Figure 4. Geometric structure of the systems at different phases.The plots from the left to the right represent stars S1, S2, S3 and S4, respectively.The red dot represents the hot spot on the surface of the second component of star No.4.

Figure 5 .Figure 6 .
Figure5.The relationship between periods (P) and fill-out ratios (f), with a linear regression line fit (blue line) for data points of W UMa systems with periods less than 0.8 d .The black circles are the data collected from48 .The red X signs are for the present work.

Figure 7 .
Figure 7. Mass-Luminosity (upper) and Mass-Radius (lower) evolutionary tracks for the four systems.The dashed lines represent the Terminal Age Main Sequence (TAMS), while the solid lines correspond to the Zero Age Main Sequence (ZAMS).These tracks were extracted from 57 .

Figure 8 .
Figure 8. Period Angular momentum (upper) and Period Total mass (lower) plots for the four systems in comparison with 58 .

Table 1 .
Names, coordinates and color index (J-H) of the four studied systems.

Table 3 .
The average Time of minima, cataloged period, and period from the current work.

Table 4 .
Physical parameters of the four systems extracted from the light curve analysis.

Table 5 .
Absolute masses and luminosities of primary and secondary components of W UMa systems for illustrating the M-L relation in the current work.

Table 6 .
Absolute (global) parameters for each system.